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Continuous weak measurement and feedback control of a solid-state charge qubit: 
physical unravelling of non-Lindblad master equation 

Shi-Kuan Wang, Jinshuang Jin * and Xin-Qi Li^ 
State Key Laboratory for Superlattices and Micro structures, Institute of Semiconductors, 
Chinese Academy of Sciences, P.O. Box 912, Beijing 100083, China 
(Dated: February 6, 2008) 

Conventional quantum trajectory theory developed in quantum optics is largely based on the 
physical unravelling of Lindbald-type master equation, which constitutes the theoretical basis of 
continuous quantum measurement and feedback control. In this work, in the context of continuous 
quantum measurement and feedback control of a solid-state charge qubit, we present a physical 
unravelling scheme of non-Lindblad type master equation. Self-consistency and numerical efficiency 
are well demonstrated. In particular, the control effect is manifested in the detector noise spectrum, 
and the effect of measurement voltage is discussed. 



I. INTRODUCTION 

Quantum trajectory theory has been developed and in- 
tensively applied in quantum optics 0, 0, H, Q, H, H, 0] • 
However, the theory is largely based on unravelling of 
Lindblad-type master equation, which has clear physical 
interpretation Q. For non-Lindblad type master equa- 
tion |8|, or even the non-Markovian dissipative dynam- 
ics |9l. I lOt Hi! [l2| , unravelling in terms of stochastic dif- 
ferential equation has also been established for various 
purposes. Nevertheless, these schemes lack physical in- 
terpretation, having thus only mathematical meaning. 

In solid-state system, the quantum measurement (T3 , 
Il4l [l~5L [l6l ll~7| of solid-state qubit and feedback control 
[18l Il9 . l20j have been an extensively studied subject in 
recent years, being largely stimulated by the prospect of 
solid-state quantum computing. The theoretical descrip- 
tion of this solid-state qubit measurement problem was 
developed originally in terms of the "n" -resolved mas- 
ter equation [lj, [l4| > where "n" stands for the number 
of electrons passed through the measurement apparatus. 
Alternatively, Bayesian formalism (la and also the con- 
ventional quantum trajectory theory [171 ] were developed 
for this solid-state measurement setup. Being of inter- 
est, it can be shown that all these three approaches are 
precisely equivalent to each other . 

More specifically, for the setup of a pair of coupled 
quantum dots (CQDs) qubit measured by a quantum 
point contact (QPC) detector, which has become an ex- 
perimentally studied system [21[ , we found that the above 
mentioned theories were restricted to the limit of large 
measurement voltage across the QPC [22|, [23|. At finite 
voltages, i.e., as the measurement voltage is comparable 
to or not much higher than the qubit 's intrinsic energy 
scale, the measurement dynamics is governed by a non- 
Lindblad type master equation [22I ]. In this regime, the 
"n" -resolved master equation was also developed to study 



the readout characteristics [23|]. In the present work, we 
further extend it to a quantum trajectory theory, which 
conditions the state evolution on the entire measurement 
records. It is well known that this kind of description 
is essential to the quantum feedback control and other 
possible applications. 

The paper is organized as follows. In Sec. II we first 
present a brief description for the setup under study, then 
outline the "n" -resolved master equation for the mea- 
surement. In connection with the conditional state evo- 
lution, physical unravelling scheme is constructed and 
practical Monte Carlo simulation is carried out. Also, 
for completeness and latter use, a generalized quantum- 
jump approach is formulated for the calculation of noise 
spectrum. For the sake of brevity, some complicated ex- 
pressions and mathematical details are put in two Ap- 
pendices. In Sec. Ill we apply the developed formalism 
to the study of feedback control, based on a suboptimal 
feedback algorithm. Finally, conclusion is presented in 
Sec. IV. 



II. CONTINUOUS MEASUREMENT OF A 
CHARGE QUBIT 

A. Model Description 

Let us consider a solid-state charge qubit measured 
by a nearby mesoscopic detector. In this work, as well 
studied in literature [H, EE EE E3 > the charge qubit is 
modelled by a pair of coherently coupled quantum dots 
with an extra electron in it, and the mesoscopic detector 
can be a mesoscopic quantum point contact (QPC). The 
Hamiltonian of this qubit-plus-detector setup reads 
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In this decomposition, the free part of the total Hamilto- 
nian, Hq, contains Hamiltonians of the measured qubit 
Hqb and the QPC's reservoirs. The operators ct(ck) 
and d k (dk) are, respectively, the electronic creation (an- 
nihilation) operators of the QPC's left and right reser- 
voirs. The qubit states \a) and \b) correspond to the 
electron locating in the left and right dots. Introducing 
e = (e a — e&)/2 and taking (e a + e b )/2 as the reference en- 
ergy, the qubit eigenenergies read E\ = ^/e 2 + d 2 = A/2, 
and E% = — Ve 2 + d 2 = —A/2. Correspondingly, the 
eigenstates are |1) = cos||a) + sin||6) for the excited 
state, and |0) = sin||a) — cos for the ground state, 
where 9 is defined by cos 8 — 2e/A, and sin 9 = 2d /A. 

For the sake of simplicity, we assume T kq = T and 
Xkq = X: i- e -j the tunneling amplitudes are real and 
reservoir-state independent. Corresponding to the qubit 
states | a) and \b), the stationary detector currents read 
I a = 2ng L g B {T + x) 2 V, and I b = 2ng L g R T 2 V, respec- 
tively, where gL(R) is the density of state of the left (right) 
reservoir, and V is the measurement voltage. Physically, 
AI = I a — If, characterizes the detector's response to 
the qubit electron's location in the CQDs. The detec- 
tor is said in the so-called weakly responding regime if 
AI <C Iq = (I a + h)/ 2. In this work we assume this 
regime, which enables us to ignore the individual electron 
tunnelling events and treat the current as a continuous 
diffusive variable. 



B. "n" -Resolved Master Equation 

The so-called "n" -resolved master equation is obtained 
by partially tracing out the detector's microscopic de- 
grees of freedom but keeping track of the number "n" of 
electrons that have tunnelled through the detector dur- 
ing the time period (0,i) [H, 0, Hia]. Originally, it was 
derived in Ref. [l3| from the many-body Schrodinger equa- 
tion at zero temperature and in the large measurement 
voltage limit. Later, it was proved that this approach 
is completely equivalent to the Bayessian approach [HI 
and the quantum trajectory theory [l7| . which share the 
same Lindblad type master equation and its unravelling. 

Alternatively, under arbitrary (i.e. not high enough) 
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For brevity, the explicit expressions of the coefficients 
a i(2), &i(2) and Ci(2,3) are ignored here and are put alter- 
natively in Appendix A. 

Formally, we reexpress the Fourier-transformed master 



voltages, it was found that this measurement problem 
cannot be described by a Lindblad type master equation 
f22[ l23l| . In particular, a non-Lindblad type master equa- 
tion was derived [22j], and its "n" -resolved counterpart 
reads 

-g( + ) p (" +1 )g + H.c.}. (2) 

Here, Q = T + X \a)(a\, Q - Q<+» + Q<~\ Q (±) = 
&W(C)Q, and C , ( ± )(£) = /^^(ije-^. C^\t) 
are the reservoir electron correlation functions. Under 
wide-band approximation for the QPC reservoirs, the 
spectral function C^(C) can be explicitly carried out 
as C(±)(£) = 2k 9l9r [x/(1 - e- x/T )] x= _ CTV , where T 
is the reservoir temperature (in this work we use the unit 
system of Ti = e = ks = 1). It is of interest to note that 
the Liouvillian operator C in C^(C) contains the infor- 
mation of energy exchange between the detector and the 
qubit, which correlates the energy relaxation of the mea- 
sured qubit with the inelastic electron tunnelling in the 
detector. Note also that in the derivation of the above 
"n" -resolved master equation we did not make assump- 
tion of large bias voltage across the QPC detector. At 
large voltage limit, i.e., the bias voltage is much larger 
than the internal energy scale of the qubit, the spectral 
function CK±)(£) ~ C^(0), and Eq. © reduces to the 
result obtained in Ref. [l3|. 

Formally, we rewrite Eq. © as 

P {n) = -iCp {n) - Kp {n) + + TZ 2 p (n+1) , (3) 

where 1Z, IZi and IZ2 are superoperators defined in accord 
with Eq. ©. To solve this infinite number of coupled 
equations, we perform the discrete Fourier transforma- 
tion, p(k,t) = £„ em V n) (0, yielding 

p(k,t)^[-iC-TZ + e lk n 1 +e- lk n 2 ]p(k,t). (4) 

Explicitly, in the localized dot-state representation 
{\a), \b)} we obtain 



(5) 



I 

equation as p(k,t) = M.(k)p(k,t) 1 and the solution reads 
p(k,t) = e A1 (' £ ^ t ~ to )p(fc, t ). Note that we are concerned 
with the "n" -resolved state evolution from to to t, i.e., 
the counting of "n" starts from the moment to. We thus 
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have p (n \t ) = p(t Q )S n fl, and p(k,t ) = p(t ). With the 
knowledge of p(k,t), the inverse Fourier transform gives 

pM(t) = J2^ mk p{k,t) = Y / e~ mk e M(k)(t - to) p{t () ).(6) 

k k 

Strikingly, we can introduce a propagator for the state 
evolution, U{n,t) — e - mfc e - A/( ( fe )*. Since this propa- 
gator is completely determined by the dynamic structure 
of the master equation but does not depend on the ini- 
tial state, we can numerically evaluate it by a "one-time 
task" via such as the fast Fourier transformation. This 
feature leads to a very efficient Monte Carlo simulation 
for the measurement-history conditioned evolution (i.e. 
the quantum trajectory simulation). 

C. Monte Carlo Simulation for Conditional 
Evolution 

It seems that p( n \t) contains less information about 
the measurement record (history) than the conditional 
state p r jt) in the conventional quantum trajectory the- 
ory [l7| , since it only implies that totally there have been 
n electrons passed through the QPC junction during the 
specified time period. However, if we make successive 
readout for the electron numbers "rife" passed through 
the detector during the time interval (tk-i,tk), we ac- 
tually record the measurement current I c (t) of a single 
realization. After each time of reading out "rife" , the 
statistically mixed state p^ n '(tk) with any possible "n" 
would "collapse" to a normalized state p^ nk \tk) with def- 
inite "n — rife". The set of records {rife : k = 1, 2, • • • } 
corresponds to the current I c (t), and the set of states 
{p( nk ^ (ifc) : k = 1, 2, • • • } is nothing but the conditional 
state p c {t) in the quantum trajectory theory. 

This state-update procedure based on the "n" -resolved 
master equation was introduced in Ref. ITEl where its ex- 
act equivalence to the Bayessian and quantum trajec- 
tory theories was analytically proved. That is, the condi- 
tional master equation can be re-derived based on the 
"n" -resolved master equation together with the above 
"collapse" idea. However, we found that this can be done 
only in the large voltage limit which leads to a Lindblad- 
type master equation [TBI, [T3] - For arbitrary voltage, 
rather than deriving a conditional master equation to 
describe the measurement-record conditioned evolution, 
we would like here to develop an efficient numerical un- 
ravalling scheme which has the advantage of being appli- 
cable to non-Lindblad type master equation as studied 
in this work. 

More quantitatively, let us consider the state evolution 
during [tj,tj + r]. That is, starting with a definite state 
at tj, say p(tj), the state p^ nj \tj + r) at tj + r can be 
calculated via 

p^\t 3 +T)=U{n 3 ,r)p{t 3 ). (7) 
If the measurement is made but the result is ignored, the 



(mixture) state is described by 

rij rij 

where Pr(rij) = Tr[p( n ^ (tj+r)] stands for the probability 
having rij electrons tunnelled through the detector, and 
p c (nj,tj + t) = p( nj \tj + r)/Pr(rij) is the normalized 
state conditioned by the definite number of rij electrons 
observed passed through the detector. 

The second equality of Eq. © implies that if we 
stochastically generate n 3 according to the probability 
Pr (rij) for each time interval [tj,tj + r], step by step 
from to to t, and "collapse" the state definitely onto 
p c (nj,tj + t), i.e., p c (tj + r) = p^fa + T )/ Pr (^i) = 
U (nj,T)p(tj)/Ti[U (rij, r)p(tj)], we have in fact simulated 
a particular realization for the selective state evolution 
conditioned on the (continuous) specific measurement re- 
sult. The simple ensemble average over a large num- 
ber of particular realizations of p c (t) recovers the un- 
conditional state p(t). Obviously, this unravelling scheme 
is completely equivalent to the spirit of the conventional 
quantum trajectory theory, despite that in this context 
we are unable to derive an explicit stochastic differential 
equation to unravel the underlying non-Lindblad master 
equation. However, they have precisely the same physical 
meaning. 

To stochastically generate nj according to the prob- 
ability Pr(rij), two procedures are adopted as follows: 
(i) Based on the ensemble average current I(t) — 
Re{Ti[Q p(t)Q]} , which was derived in Ref. [H, the out- 
put current in particular measurement realization reads 

I c (t) = Re{Tr[g Pc (t)Q]} + £(t). (9) 

The first term in this equation is related to the qubit 
dynamics of conditional evolution. The second noisy 
term £(t) originates from the detector's intrinsic noise, 
which is a Poisonian variable in the regime of point pro- 
cess, and a Gaussian variable in the diffusive regime. In 
the latter case, has zero mean value, and the spec- 
tral density — 2/oCoth^. At zero temperature and 
large voltage limit, this treatment recovers the existing 
result of q uantum trajectory theory and Bayesian ap- 
proach [HI, [l7|], i-e., ijt) = Pc,aa(t)Ia + Pc,bb(t)h + £(*), 

with S% = 2Iq. (ii) Straightforwardly, in our simula- 
tion we relate the stochastic electron number rij with 
I c (t) via nj = $l s+T dt'I c (t') = I c (tj)T + dW(tj), where 

I c (tj) = Re{Tr[Qp c (tj)Q]}, and dW(t 3 ) is the Wiener 
increment during [tj , tj + r] . 

In Fig. 1 we plot a comparison of the ensemble average 
of the Monte Carlo simulation (over 500 quantum trajec- 
tories) with the result directly given by the unconditional 
master equation. The excellent agreement shows the va- 
lidity and efficiency of the proposed unravelling scheme. 
In this context, two points are likely to be highlighted: (i) 
the measurement voltage considered here is moderately 
finite, but not the high voltage limit [IE G3; (ii) tnc 
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FIG. 1: Ensemble average of the Monte Carlo simulation over 
500 quantum trajectories (dash-dotted line) versus the result 
directly given by the unconditional master equation (solid 
line) . It is assumed that the initial state of the qubit \ip) = a) . 
The relevant parameters are: T = 20Q, \ = 0.7f2, V = 0.5S1, 
e = 0.25O, T = l.Ofi, and g L(R) = l/V2nSl. 



corresponding non-Lindblad master equation is unrav- 
elled physically, having the same physical interpretation 
as provided by Wiseman et al @ . 

In Fig. 2 we show the main features of the conditional 
state evolution, (i) Assuming that the coherent coupling 
is switched off (fi = 0), in Fig. 2(a) we illustrate the wave- 
function collapse of a pure state = l/\/2(\a) + \b)) 
under measurement. This feature has deep implication 
in understanding the measurement postulate in quantum 
mechanics, which has been highlighted by the concept of 
gradual collapse for a typical solid-state two-level state 
under (weak) measurement, as discussed in particular by 
Leggett [2~ij ]. The reason for this gradual collapse is the 
weak coupling and the finite noise of the detector, which 
make the quantum measurement need some time until 
acceptable signal-to-noise ratio is reached, (ii) On the 
other hand, if ^ 0, the ideal measurement will lead 
to the gradual purification of the qubit state starting, for 
instance, with a completely mixed state. This is shown 
by the revival of the coherent Rabi oscillation in Fig. 
2(b). (iii) As shown in Fig. 3(c), with the increase of 
the measurement strength (|^|/0), the duration time on 
each qubit state is enhanced, while the switching time 
between them is reduced. This is an obvious signature 
of the quantum Zeno effect, appearing in the regime of 
gradual but not the conventional instantaneous collapse, 
(iv) In the conditional dynamics, the output current I c {t) 
can basically follow the conditional qubit state, as shown 
by Fig. 2(d). Due to the intrinsic noise of the detector, 
here the filtered current, i.e., I c (t) = J^ +At I c (t )dt , 
is plotted, where At is the "filtering window" . 



FIG. 2: The main characteristics of the conditional state 
evolution under continuous measurement: (a) gradual local- 
ization (fi = 0) from an initial superposition state — 
1/V2(|a) + 1 6}); (b) gradual purification from a completely 
mixed state; (c) zeno effect under relatively strong continu- 
ous measurement; and (d) conditional state evolution (thick 
line) versus the filtered output current (thin line), which 
is obtained by J ' I c {t)dt , with the filtering window 
At = 0.2/Q,. The parameters (T, x, V, e, T) are adopted as 
(in unit of fi): (a) (20.0, 0.13, 3.0, 0.25, 1.0), (b)(20.0, 0.7, 3.0, 
0.25, 1.0), (c)(25.0, 4.0, 3.0, 0.25, 1.0), and (d)(25.0, 4.0, 3.0, 
0.25, 1.0). The density of states gi(m = 1/V2tvQ. 



D. Noise Spectrum 

In the same spirit of conventional quantum trajectory 
theory 0, [T3|, the present unravelling scheme also pro- 
vides a natural way to calculate the output power spec- 
trum. The details of derivation is referred to Appendix 
B, here we simply present the resultant expression of the 
output current correlator, which reads 



Kj{t) = E[I c (t + r)I c (t)}- E[I c (t- 
= Tr[Ue CT Up{t)} - Tr[Up(t - 
+Tr[U 7 p(t)}5(r), 



-r)]E%(t)] 
T)]Tr[Up(t)] 



(10) 



where U = £\ n^(n, tft)/dt, W = ^ n n 2 W(n,rft)/dt, 
p(t + t) = e p(t). For stationary state Kj(t) = 
Tr[Ue CT Up(oo)} - Tr[Z7p(oo)] 2 + Tr[W 7 i o(oo)](5(r). 

In practice, particularly in the presence of quantum 
feedback, the ensemble averaged evolution represented 
by e Cr can be implemented by numerically averaging the 
stochastic trajectories. In the absence of quantum feed- 
back, the above quantum trajectory approach can pre- 
cisely recover the analytic result of stationary state noise 
spectrum, obtained by using the MacDonald formula [23[ . 



III. QUANTUM FEEDBACK CONTROL 

Quantum feedback control is one of the typical means 
of quantum coherence control. In quantum optics, the 
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study of quantum feedback control has been going on 
for more than a decade [25[. However, it is a relatively 
new subject in solid states [H, EH, Ho]- In particular, 
the conditional state evolution under continuous weak 
measurement has been experimentally demonstrated in 
solid-state qubit very recently [2(| . This may pave a way 
to the quantum feedback control in solid-states. For the 
solid-state setup under present study, we now consider 
the feedback control of the qubit coherent evolution, by 
unravelling the underlying measurement dynamics that 
is in general governed by non-Lindblad master equation. 

The basic idea is to convert the measurement infor- 
mation in the output current I c (t) into the evolution 
of qubit state p c {t). By comparing p c (t) with the de- 
sired state pd(t), their difference is then employed to 
modify the qubit Hamiltonian in order to reduce their 
difference in next step evolution. Specifically we con- 
sider a symmetric qubit (i.e. e = 0). The desired state 
is \4>d{t)) = cos(f2i)|a) — ism(Q£)\b). In real-time feed- 
back control, each successive feedback acts only for an 
infinitesimal time interval At. In the so-called Bayessian 
state-estimate-based feedback, the suboptimal algorithm 
is desirable [27ll28|. That is, the algorithm is constructed 
such that the state evolution in each infinitesimal time 
step will maximize the fidelity of the estimated state with 
the desired (target) state. In more detail, as far as the 
term related to the feedback Hamiltonian is concerned, 
the state p c (t + At) is given by 

Pc (t + At) = p c (i)-i[H fb ,p c (t)]At 

-^[H fb7 [H fb , Pc (t)}}(At) 2 + -- - .(11) 
The fidelity of this state with the target state reads 

(Mt)\pc(t + At)\Mt)) 

= (Mt)\Pc(t)\Mt)) - i(Mt)\[H fb, Pc {t)]\4>d(t))At 

-^(Mt)\[H fb AH fb ,Pc(t)}}\Mt))m 2 + -- - ■ (12) 

To optimize the fidelity, one should maximize the coef- 
ficient of At, which is the dominant term. Similar to 
other control theories, the maximization must be sub- 
ject to certain constraints, e.g., the restriction on the 
maximum eigenvalue of Hf b , the sum of the norms of 
the eigenvalues, or the sum of the squares of the eigen- 
values, etc. Physically, these constraints stem from the 
limitation of the feedback strength or finite Hamiltonian 
resources. Here we adopt the last type of constraint, 
namely, Tr[Hj b ] < p . Under this constraint, the feed- 
back Hamiltonian can be constructed in terms of 

H fb = i\[\Mt))(Mt)\,Pc(t)] , (13) 

where A = y^^y, with a = {ip d (t)\p c (t) 2 \ip d (t)) , and 

b = (Mt)\pc(t)\Mt))- 

Combining the above feedback Hamiltonian with the 
previously developed state unravelling scheme, the es- 
timated state p c (t) can be straightforwardly calculated, 
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FIG. 3: Feedback control of the coherent oscillation of qubit 
state, resulting from feedback strengths A = 0.5 (thin solid 
line), 1.0 (thick solid line), and 3.5 ( dashed line). The pa- 
rameters are: T = 20.0O, % = 0.7fi, V = 3.0fi, e = O.Ofi, 
T = l.Ofi, and g L(R) = l/y/2nSl. 



leading to the state propagation in the presence of feed- 
back. Figure 3 shows the control result with quantum 
feedback, where the ensemble-average has been made 
over large number of Monte-Carlo simulated trajecto- 
ries. Here the measurement voltage is quite moder- 
ate, i.e., V — 30, which is beyond the theoretical 
description in larg e- voltage limit as previously studied 
[II EE E3, 0, EII E3|- We observe that the control ef- 
fect is evident: by increasing the feedback strength A the 
measurement induced back-action can be largely elimi- 
nated, and the desired coherent oscillation of the qubit 
can be maintained for arbitrarily long time. 

It will be of interest to compare the quantum measure- 
ment in the presence of feedback to the well-known quan- 
tum non-demolition (QND) measurement. For the solid- 
state qubit, elegant schemes of QND measurement have 
been proposed very recently [2{J,|3(|. Here, for the qubit 
measurement under consideration, we demonstrate that 
the QND measurement is equivalent to the usual back- 
action-present measurement plus quantum feedback. In 
Fig. 4 we present the calculated output power spectrum 
of the QPC detector. The peak at u = Q indicates the 
coherent oscillation of the qubit. In the absence of feed- 
back, it has been shown that the peak-to-background ra- 
tio cannot be larger than 4 [l6j] , due to the back-action of 
measurement. In the presence of feedback, however, we 
obtain very sharp peak here which indicates almost ideal 
coherent oscillations. Theoretically, since no steady-state 
is available in the presence of feedback, the start time of 
the qubit evolution is chosen as the initial time of the 
current correlation function, and the noise spectrum is 
the Fourier transform of the correlate function with re- 
spect to the later time (difference). Experimentally this 
feedback-induced sharper peak can be employed as an in- 
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FIG. 4: Narrowed coherence peak in the noise spectrum as an 
indicator for the feedback effect (Inset: the coherently oscillat- 
ing correlation function). The parameters are : T = 10.0O, 
X = 0.2fi, V = 3.00, e = 0.00, T = 0.5O, A = 15, and 
9l(r) = 1/V2ir0. 



FIG. 5: Effect of measurement voltage on the feedback con- 
trol. Here the quantity synchronization degree, which is de- 
fined as D = 2{Tr[p c (t)pd(t)]) — 1, with (• • • } meaning the 
average over time, is employed to characterize the control 
quality. Note that for perfect control the synchronization de- 
gree is unity (i.e., D — 1). The parameters are: T = 20.0O, 
X = 0.7O, e = O.Ofi, T = LOO, and g L(R) = l/V2nQ. 



dicator for the feedback effect in practice. We expect that 
this kind of experiment can be performed in the not-far 
future. 

Finally, we address the effect of measurement voltage. 
Figure 5 shows the synchronization degree of the feed- 
back versus the measurement voltage. It is of interest to 
notice that there exists an optimal measurement voltage 
for relatively small feedback strength. This is because 
for larger voltage the back-action is relatively too strong, 
while for smaller voltage the information of the measured 
state cannot be extracted out efficiently. As a result, the 
turnover behavior of the synchronization degree versus 
voltage is found. However, with the increase of the feed- 
back strength, the strong back-action can be eliminated 
more efficiently by the feedback. In this case, the syn- 
chronization degree does not decrease considerably as the 
measurement voltage increases, as shown in Fig. 5. 



IV. CONCLUDING REMARKS 

To summarize, in the context of continuous quan- 
tum measurement and coherence control of solid-state 
charge qubit, we have presented a unravelling scheme 
for the non-Lindblad type master equation. Based on 
it, we also constructed an efficient method to calculate 
the noise spectrum, which can be regarded as a general- 
ization of the standard quantum jump theory developed 
in quantum optics. Despite the absence of analytic for- 
malism, the numerical implementation in practice was 
demonstrated to be straightforward and efficient. Illus- 
trative application was further contributed to the quan- 
tum feedback control under arbitrary measurement volt- 



ages. The detector noise spectrum under feedback was 
calculated, and its narrowing clearly reflected the con- 
trol effect. Also, the effect of measurement voltage was 
discussed. 

The present study has been focused on the setup of 
double-dot qubit measured by quantum point contact. 
However, for other solid-state setup such as charge qubit 
measured by single-electron-transistor, similar unravel- 
ling scheme can be constructed. Owing to the fact that in 
general (e.g. in the presence of many-body Coulomb cor- 
relations) the measurement dynamics is not governed by 
Lindblad type master equation, the present "n" -resolved 
master equation based unravelling scheme seems quite 
desirable. Finally, in spite of the various unravelling 
schemes for non-Lindblad type master equation or even 
for non-Markovian dissipative systems, to our knowl- 
edge all of them are largely mathematical. Therefore, 
the present physical unravelling scheme is of interest and 
valuable, which may find applications in the field of solid- 
state quantum information. 
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APPENDIX A: MATRIX ELEMENTS OF M{k) 
The coefficients in Eq.(5) read 



Oi = 


Qaa [Qaa 


- e ik Q& - e~ ik Q£>], 


(Ala) 


«2 = 


Qaa r/\ 
Y~ iQ ab 




(Alb) 


h = 


— Qbb{Qbb - 


JkA(-) P - ik n( + h 


(Ale) 


b 2 = 


Qbb r^j 

2~ [Qba 




(Aid) 


Cl = 


7^ [(QaaQt 


jo + QbbQbb) - e lk (QaaQ~l b ^ 


+ QbbQ ( J) 



C2 = 
C3 = 

where 
rj(±) 



Qaa 

and 



-e^iQaaQ^+QbbQiV)}, 

i, 



(Ale) 



2 [QaaQab - e 4fe Q b6 QL } - e-^Q^] , (Alf ) 

~[Q66Q6a - ^QaaQ^ - e^Q,.^], (Alg) 



= C< ± H0)[T+ 1 X (1 + cos 2 0)] + X A± sin 2 0, 
- C( ± )(O)[T+i X (l + sin 2 0)]- X A ± sin 2 ^ 
< 3ib ) = ^xC ,(±) (O)sin6»cos6» + xsin0(A ± - A±cos6»), 
\xC {±) (0) sin 6 cos - x sin <9(A± + A± cos 9), 



T+ Xi Qbb — T, Q ab — Q b 



A± = J[(7( ± )(-A) + C( ± )(A)], 
A± = ^'(-AJ-C'^A)]. 



APPENDIX B: NOISE SPECTRUM 



In this appendix, along the line of the conven- 
tional quantum jump theory, we extend the method 
of noise spectrum calculation to the unravelling ap- 
proach developed in this work. Consider the corre- 
lation function E[dn(t + r)dn(t)], where E[dn(t)] = 
J2 ni niTr[U(m,dt)p(t)]. First, for the case r > 0, 



E[dn(t + r)dn(i)] 

= ^niProb[dn(t) = ni]E[dn c (t+T)\ dn ^ 



(Bl) 



where Prob[dra(t) = m] = Tr[W(m, dt)p{t)}. At time 
t + dt a definite value n\ is picked out and the qubit 
state undergoes an immediate collapse, i.e. p c (t + dt) = 
U(ni,dt)p(t)/Prob[dn(t) = m]. During the time period 
[t + dt, t + t], one can ignore the measurement records 
owing to the ensemble nature of the correlation function, 
therefore 



E[dn c (t + r)| dn ( t)=Tll ] 

= Zr»T±MT»,dt)^-*) J%^ . ]. (B2) 
^ Prob[dn(i) = mj 



To the leading order of dt we have 



E[dn(t + T)dn{t)] = J2 n 2 TT[U(n 2 ,dt)e CT x J^L^^_ jmTr^K, eft)p(t)] = Tr[We £ ^p(t)]dt 2 . (B3) 



Here we have introduced W = ^2 n nL((n, dt)/dt. Next for 
t = we have 

E[dn(t) 2 } = ^n 2 Prob[dn(i) = n] 

n 

= ^VTr[W(n,dt)p(t)] 

n 

= Tr$Fp(t)]dt, (B4) 

where W — ^ n n 2 U(n,dt) / dt. For short time 
r this equal-time correlation will be dominant, and 
E[ dn(t d + r) - ^1] can be treated as ^-correlated noise for 



a suitably defined 5 function. We thus obtain 

K/(t) = E [ dn{t + T) 1 E[ dH{t + T) ]^[ rfn(t) 1 
di dt dt 

= Tr[Z7e £T Z7p(i)] - Tr[Z7/9(t + r)]Tr[Z7p(i)] 
+Tr[W 7 p(t)]<5(T). (B5) 

Finally, it should be noted that [dn{t)\ 2 ^ dn(t) in 
our above treatment. This differs from the conven- 
tional quantum jump theory where the stochastic number 
dn(t) = or 1 in the point process. 
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